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ABSTRACT 

The navigational requirements of Galileo as it swings by Mars 
[flyby distance is 275 km from the Martian surface; 25 km (lo) ] are 
going to be met with interferometric angular measurements (VLBI) and 
range and range-rate measurements. Like VLBI, dual spacecraft differ- 
enced range is less sensitive to Mars ephemeris errors and tracking 
station location errors than conventional range and Doppler. Similarly, 
differenced range provides angular information about the separation 
between the Mars Viking Lander I and the Galileo spacecraft. In covariance 
studies, dual spacecraft range coupled with conventional range and Doppler 
is shown to estimate the Galileo-Mars flyby distance to better than 10 km 
(10) which is comparable to the VLBI performance. For the Galileo-Mars 
flyby, dual spacecraft differenced range promises to be an excellent 
backup to VLBI if the Mars Viking Lander remains operational. 


This paper presents the results of one phase of 
Jet Propulsion Laboratory, California Institute 
NAS7-100, sponsored by the National Aeronautics 


research carried out at the 
of Technology, under Contract 
and Space Administration. 
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I. INTRODUCTION 

Galileo, a NASA spacecraft to be launched in 1984 by the Space Shuttle/IUS 
launch vehicle, will travel from earth to a Mars flyby on an ultra fast trajec- 
tory: a flight time of less than 100 days. On such a high acceleration tra- 
jectory (Fig. 1), conventional Galileo radiometric tracking data, 2-way range 
and Doppler, can establish the heliocentric position of the probe to a standard 
deviation of II km. 

The Mars ephemeris has an additional 40 km * in- track* position uncertainty 
such that the Mars-Galileo relative position uncertainty exceeds 4O km (la). 

It is essential to know the Mars— Galileo relative position to better than 25 km 
(1<7) • The closer Galileo can be flown past Mars, the smaller the Galileo rocket 
maneuver that will be required to -send Galileo on to the Jupiter system (Fig. 2). 
The AV requirement of Galileo's rockets increases 40 m/s per 100 km increase in 
the flyby distance (Fig. 3). The Mars flyby is being used to provide a con- 
trolled acceleration to the Galileo spacecraft. 

Deep space probes, such as Geillleo, are tracked and navigated from earth. 
That is, a radio carrier is beamed to a distant space probe. The probe trans- 
ponds the radio tone back to earth. The frequency difference between the 
earth transmitted and received signal is the Doppler shift - a measure of 
the spacecraft radial velocity. Modulation placed on the radio carrier is 
used to measure the light time separation between earth transmission to 
and reception from the spacecraft. These conventional radio metric data 
types, Doppler and range, measure in the radial direction only. 
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Because of topocentric parallax, at any instant of time the radial velocity 
of a probe is different at each terrestrial tracking station. This uniqueness 
of Galileo and Viking'^ range and range-rate, that is dependent on topocentric 
position, permits the determination of the relative Galileo-Viking angular 
separation. 

In covariance studies, Galileo conventional range and Doppler could 
estimate the Galileo-Mars flyby distance to 40 km (a) while a combination 
of doubly differenced range which exploits the relative topocentric parallax 

in conventional range and Doppler yields a standard deviation of 
less than 10 km and does so 25 days before Mars encounter (Fig. 4) . 

Galileo Project plans call for the Galileo spacecraft to flyby Mars 275 km 
(o^ = 25 km) above the planet’s surface (Ref. 1). To achieve this accurate 
flyby two new technological advances must be accomplished: one, the Mars 

ephemeris must be improved to better than 25 km (a) and this effort is in pro- 
gress; two, a wide-band Very Long Base Interferometry technology must be deve- 
loped that will permit the Galileo spacecraft and Mars trajectories to be 
defined in a quasar inertial reference frame. This latter effort is underway 
also and offers not only a means to reduce the Galileo-Mars relative trajectory 
errors but VLBI cancels the preponderance of the Beep Space Station (DSS) 
Iccatica effects on orbit determination. 


* The Viking Mars Lander I softly touched down on the Martian surface on 
July 4, 1976, and it still functions. It is expected to be operational 
in the Galileo era. 
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As circumstances are now. the current Mars Epheraeris and DSS locations 
uncertainties limit the Galileo~Mars relative navigation such that a 0^ 5 25 km 
is not achievable with Galileo radiometric range and doppler alone. Galileo- 
Viking doubly differenced range provides a promising approach to Galileo's 
navigation objectives independent of an improved Mars epheraeris or a new VLBI 
technology. It does require the survival of the Viking Lander, however, in 1984. 

II. DOUBLY DIFFERENCED RANGE DEFINITION 

Figure 5 shows two Deep Space Station (DSS) tracking first one space- 
craft and then the other. Thus, four range measurements are obtained and 
although the order of the range measurements taken in Figure 5 are DSS-1 to 
Viking, DSS-1 to Galileo, DSS-2 to Viking and DSS-2 to Galileo, the order ts 
arbitrary. 

With a restricted view to a single spacecraft, it is easy to show that 
the relative topocentric range (Fig. 6) involving 2 DSS is 
Ap ** AZ sin 6 + AL cos 6 

where 

Ap « p2 " Pj^ 

AZ = north-south projection of the DSS baseline on that plane 
possessing the baseline and the spacecraft 

AL * east-west projection of the baseline 

or 

AL ® AX cos (a - LST) 
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with 


AX being a linear separation between the DSS in the earth 
equatorial plane 

LST *= local sidereal time at the X of the baseline 


( *DSSj ) 


/2 


Now If the relative range, Ap, from 2 spacecraft are combined in a second 
difference 


A p - APg - Ap^ 


G signifies Galileo 
V signifies Viking 


A^p - AZ 1 

cos 

«G “1 



- AX 

sin 

ag - LST) 

cos 6_ Aa 

Car 

(1) 

+ 

cos 

(“g ‘ 

sin 6 A5 


is a function 

of 

the relative 

plane-of-sky coordinates of the 

two 


spacecraft and the baseline projection onto the plane-of-sky . It*s sensitivity 
to the Mars ephemeris is less than that of Ap or p. 

8Ap, 


^a2 

3A p 


3Ap^ 


^^state) ^c^(state) ^<^(state) 


but 


3A^p ^ 3Ap 

9 (Galileo State) * 3(Galileo State) 

2 

Specifically, A p is 20% (2 months before Mars encounter) to 50% (at encounter) 

less sensitive to the Mars ephemeris error than Ap* as is shown by the RSS of 

u 

2 

Ap and A p partials with respect to the heliocentric position of Mars (Fig. 7) . 
In Figure 7 there are three graphs, one for each baseline used in the study. 
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DSS43 is located in Woomera, Australia. DSS63 is located in Madrid, Spain, 
and DSS14 is at Goldstone, California. The DSS approximate spherical coor- 
dinates are tabulated in Table I. 

Since the DSS are separated in longitude, from 94° to 154°, Mars is in 
view over each baseline at different times. The Viking Lander can only be 
ranged in the cool morning Martian hours and only ranged once per Martian day. 
Thus, as indicated in Figure 7, the Viking Lander can be ranged about 50% of 
the days that Galileo is in flight. Table II presents the 43 different occa- 
sions. Each baseline can range the Lander for 8 to 12 days repetitively. 

Each baseline’s performance is not only time dependent, but is also 
governed by the alignment of the baseline with respect to the Galileo— Mars 
angular separation at encounter (Eq 1). In essence, the Viking and the 
Galileo Ap measurements provide information as to the direction of each space- 
craft with respect to the baseline but only in the direction of the baseline. 

Orthogonal to the baseline there is no information. And, of course, when Ap 

2 2 

measurements are differenced to obtain A p, A p defines the earth centered 
angular separation between the two spacecraft only in the baseline direction. 
Figure 8 shows the baseline orientations relative to the Mars— Galileo direction 
at encounter. The DSS43 - DSS63 baseline which is approximately 4® offset, 
yields the strongest information concerning the flyby distance while the DSS63 
- DSS14(^12® offset) and the DSS14 - DSS43 ( 60° offset^ baselines provide 
progressively less information. 

2 

Table III Itemizes the theoretical error assessments of A p resulting 
from instrumentation and transmission media** From Table III it is apparent 
Philip Callahan, Jet Propulsion Laboratory, private communication 
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that it is thermal noise [galactic backgrouiid (6*K) , receiver front end elec- 
tronics (6“-ll® K) , antenna cable (3* K) , transmission media (10® K) , etc.] 

2 

when subjected to high gain that dominates the A p error budget. A full 90% 

2 2 
of the A p RSS noise is from this source. The result is that A p should 

have an RMS error of'V'2 m and 2m is the a priori standard diviation used in 
covariance study. Most of the systematic errors due to solar plasma, tropos- 
phere, ionosphere, DSS clock errors, spacecraft and station delays cancel. In 

addition, since cracking stations are used redundantly to crack both Galileo 

2 

and Viking, DSS longitude errors tend to cancel in the formation of A p (Fig. 9). 

DSS linear tain ties in the other two coordinates are of little consequence since 

2 

their effect upon A p is from one to two orders of magnitude smaller yet. 

III. THE GALILEO-MARS FLYBY DISTANCE COVARIANCE STUDY 

The covariance analysis performed in this paper allows a maximum likelihood 
estimated with gaussian errors on the observations. The assumed observations 
include two-way coherent Doppler data from the Galileo spacecraft using the 
three Deep Space Ifetwork stations continuously, one Doppler measurement every 
one hour, one range measurement from the (toldstone station every day and the 
available doubly differenced range measuretnents as shown in the Table II. Since 
the dynamic state parameters are non-linear functions of the measurements, the 
observation equations are linearized and the results obtained are based on a 
linear estimator. VRien a standard maximum likelihood estimator is constructed, 
the computed statistics based on data noise errors, do not reflect the effect 
of model errors in the solution. Thus the statistics must be adjusted to 
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account for these effects- 

The measurement equation can be written in the form 
2 « Ax + Cp -h e 

where Z is the vector of measurements, x the vector of est Lmated parameters, p 
the vector of model parameters whose effects on the estimated parameters are 
to be investigated and e the measurement errors. A weighted least squares 
estimator of x can be obtained by (Bryson and Ho, 1969) 

X - ^ a'^P'^z 


with the assumption that p is a random vector of zero mean with covariance P^, 
E(5) * 0, Cov (5) » P and E ( pe^) “ 0 and the covariance of x is given by 


^c „ (c\-P+P A^P"^CP C^P"^AP^ 

P « Cov \ X / X X C X 

X 


where P « (a*^P ^A ) is the voir covariance matrix. The matrix P ^ is known as 

X X 

the 'consider* covariance matrix and the matrix A and C are the partial deriva- 
tives of the measurements with respect to the estimated and the consider para- 
meters. Both the Galileo orbital state and Mars ephemeris parameters are trea- 
ted as estimated parameters, and the station locations. Viking lander locations 
and Mars mass are treated as ^considered’ parameters. The apriori uncertainties 

of the parameters are given in Table IV, 

2 

In the model used to assess A p, the trajectory parameters of the Galileo 
probe was estimated in a manner that considered the uncertainties associated 
with the Mars Ephemeris, the DSS location set, the mass of Mars, the Viking Lander 
position (Table IV). 

With this parameter set and the Galileo data set (Table V) , the Galileo 
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heliocentric position can be estimated to 35 km (a). Figure 10, and this 
uncertainty steins principally from DSS location uncertainties. The Galileo 
trajectory does not sense the gravitational effect of Mars until the las t day before 
encounter. Galileo travels over a million kilometers on that last day. 

Figure 10, the ’Standard Deviation of Galileo Heliocentric Position’, 

shows a in kilometers as a function of time in days from Mars encounter or 
P 

when each simulated Galileo data arc stops. All estimates of involve data 
that starts 88 days before Mars encounter. Each estimate, following the E-85 
estimate, has an additional five days of data added to the solution. All of 
the standard deviation plots presented have this same format. 

Galileo, Mars-centered, position estimates have a standard deviation 
equal to the RSS of Galileo’s heliocentric position sigma and the. Mars ephemeris 
position standard deviation (Fig. 11a). 

Figure 11 not only exhibits the standard deviation of the Galileo flyby 

but shows the components of related to the Mars Ephemeris j Ephemeris)’ 

the DSS locations | | „oise) • 

Since Galileo is over a million kilometers away from Mars at E-1*^, Galileo 
does not see Mars gravitationally until E-2^ and any effort to utilize Galileo 
tracking data to improve the Mars ephemeris fails. Hence, the ephemeris pro- 
vides a near constant 40 plus kilometer component to o^(RSS) . 

As indicated in Figure 11, | Locations earth-probe 

distance Increases. That is, DSS angular location uncertainty in an Euclidian 
solar system results in larger and larger spacecraft linear position uncertainty 
with increased topocentric range. However, if DSS coordinates were estimated, 
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instead of considered, this procedural artifact would disappear as in Fig. 2. 

And lastly, in Figure 11, the data noise is shown to fall off with the 

square-root of the number of observations. 

2 

When these A p observations of Table II are added to the conventional 

Griileo data cf Table V, the effects of the Mars ephemeris and the DSS location 

uncertainties are reduced. Ihis should be expected since the RSS of the partials 

of the Mars position coordinates (Fig. 7) and the DSS coordinates (Fig. 9) with 
2 

respect to A p are 2 to 10 times smaller than those with respect to Ap. That is, 

2 . 2 
each A p observation is less sensitive to these error sources, but A p and Ap 

possess the same sensitivity to the Galileo-Mars relative state. Figure 12 

exhibits the ephemeris, DSS, and data noise contributions to a,. The data 

d 

2 

ensemble of A p, conventional range and Doppler yields a < 10 km (a) 25 days 

before Mars encounter. This is an improvement over conventional data reductions of 

four-fold. As can be seen in Figure 12, the correlated ephemeris and DSS loca- 

2 

txons uncertainties in each p observation cancel in the formation of A p . As 
modeled. Mars ephemeris and DSS location uncertainties still dominate the stan- 
dard deviation of the Galileo-Mars encounter distance estimate, however, their 
combined RSS contribution is less than 10 km (a) . 
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Summary 

This covariance study shows that Galileo-Mars navigation is improved 
four-fold when dual station range from both Galileo and Viking are added 
to conventional Galileo tracking data and reduced. In essence, the Mars 
ephemeris and the tracking station uncertainties are differenced out of 
the new doubly differenced range data type, to a large extent, while little 
Galileo-Mars relative state information is lost. The information content 
of doubly differenced range is analogous to that of wideband very long 
baseline interferometry and promises to be an efficient backup the Galileo 
Project plaimed VLBI. Doubly differenced range coupled with conventional 
tracking data can be used to estimate Galileo-Mars flyby distance to better 
than 10 km (o). 


Reference 

Project Galileo Navigation Requirements, ?D 625-565, JPL 19 April 1979, 
JPL Internal Document, 
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TABLE I: Tracking Station Spherical Coordinates 


DSS 

Longi tude 

Lati tude 

43 

149?0 

35?3 

63 

355.8 

-35.3 

14 

243.1 

40.3 


11J.2 






TABLE II: VIKING LANDER DIRECT LINK RANGING OPPORTUNITIES 



S. 

o 


< 


c 
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DSS 14 (Goldstone, California) 
DSS 43 (Woomera, Australia) 

DSS 63 (Madrid, Spain) 










TABLE III: Doubly Differenced NSR Error Budget 


INSTRUMENTATION: 


STATION CLOCK STABILITY (15 MIN) 

4 

CM 

STATION DELAY CALIBRATION 


-- 

SNR (THERMAL NOISE) 

200 

CM 

WAVEFORM DISTORTION 

88 

CM 

SPACECRAFT DELAY 

28 

CM 

MEDIA: 



TROPOSPHERE (25* ELEVATION) 

20 

CM 

lONOSPHE.RE (25* ELEVATION) 

6 

CM 

SOLAR WIND 

15 

CM 


RSS 222 CM 


ASSUMPTIONS: 

VIKING LANDER - GALILEO SEPARATION ~ 5® 

DATA AT OPPOSITION - 0.7 AU 

TWO STATIONS OBSERVE LANDER IN TURN APPROX. 15 MIN EACH 
SAME TWO STATIONS OBSERVE GALILEO IN TURN APPROX. 15 MIN EACH 
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TABLE IV: Gal il eo-Vi king Parameter Set 


PARAMETERS 

MODEL STATUS 

A PRIORI 

Gal il eo State 

Estimated 

0 ^ = Oy * 0 ^ * 10 ^ km; o* = ap * o* = 100 km/s 

Mars State 

Consi dered 

0 j. 1 * 10 km; 0 . * b ® 40 luff's 

radial in track 

0 . - , ^ = 70 km 

out-of-plane 

DSS Locations* 

Considered 1 

0 , = 3.0 m; 0„ = 1.5 m; 0 = 15.0 m 

A r^ r^ 

Viking Lander 
Locations 

Consi dered 

a = 10.0 m; 0 = 40.0 m; 0 = 300.0 m; 

X y z 

_3 

0 = 0 = 0=10 m/day 
X y i 

Mars GM 1 

1 

Considered 

- -J 

0 = 0.1 km'^/Sec^ 


* Tg = DSS distance from terrestrial spin-axis 
r^ * DSS distance from earth equator plane 
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TABLE V: Schedule for Conventional Data 


DATA TYPE 

(a) 

RATE 

DSS ACQUIRING 

Galileo Doppler 

1 ""/s 


1 pt/hr 

14, 43, 63 

Galileo Range 

1 km 


1 pt/pass 

14 

Start: 6 March 1984 88 days) 

Stop: 2 June 1984 20 min) 
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FIGURE 1: Solar Polar View of Galileo to Mars Trajectory 




earth at launch 



Figure 2. Heliocentric Ecliptic Pole View 
of Interplanetary Trajectory 


FIGURE 3: AV and Flyby Distance Relationship 



. S/W ‘AV 


11-19 


Minimum Altitude, km 




FIGURE 5: Range Components of Doubly Differenced Range 


Mars Lander 
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